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I.  Introduction. 

Supersonic  molecular  jet  spectroscopy  has  made  possible  the  study  of  a 
wide  variety  of  weakly  bound,  solute-solvent  van  der  Waals  vdW)  clusters  in 
the  gas  phase.  Cluster  investigations  have  enhanced  our  understanding  of 
intra-  and  intermolecular  interactions  and  potentials,  vibrational  energy 
dynamics  and  chemical  reactions,  structural  properties  of  small  aggregates  of 
solute  and'  solvent  molecules,  and  nucleation  and  growth  of  small  clusters. 
Clusters  can  also  be  considered  as  model  systems  for  condensed  phase  behavior. 
Moreover,  these  vdW  systems  can  be  thought  of  as  an  important  new  state  of 
matter  in  which  the  static  and  dynamic  properties  of  small  aggregates  of 
weakly  coupled  molecules  can  be  studied. 

vdW  clusters,  after  being  produced  in  a  supersonic  jet  expansion,  can  be 
probed  by  three  distinct  techniques:  fluorescence  excitation  (FE),  dispersed 
emission  (DE),  and  two-color  time  of  flight  mass  spectroscopy  (2-color  TOFKS). 
The  latter  technique  is  employed  most  often  in  our  studies  of  clusters  because 
it  gives  unique  cluster  identification,  brackets  the  cluster  binding  energies, 
and  elucidates  cluster  vibrational  energy  dynamics  and  vibrational 
predi3sociation. 

In  the  past  few  years,  we  have  reported  several  studies  of  vdW  clusters 


using  the  three  spectroscopic  techniques  mentioned  above. 


Cluster 


geometry,  binding  energy,  nucleation  and  growth  dynamics,  and  limits  on  the 
vibrational  energy  dynamics  and  vibrational  predissociation  times  have  been 
determined.  The  vdW  solute-solvent  clusters  investigated  initially  are  for 
the  most  part  restricted  to  aromatic  hydrocarbon  solutes  (e.f.,  benzene  and 
toluene)  and  small  alkane  solvents  (CH^,  C2H6,  and  C^Hg)  in  which  only 
one  type  of  interaction,  that  is  one  potential  form,  is  found  to  be  important 
for  the  solute-solvent  coordination.  Expanding  on  these  previous  studies  we 
are  now  exploring  solute-solvent  clusters  with  N-heterocyclic  solutes  (e.g.. 
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pyrazine  and  pyrimidine)  and  alkane  solvents  and  N-heterocvclic  solutes  with 

hydrogen  bonding  solvents  (e.g.,  water  and  ammonia).  The  initial  report  of 

this  effort  for  pyrazine  and  methane,  ethane  and  propane  clusters  has  already 
7 

appeared . 

In  this  paper  we  discuss  the  2-color  TOFKS  study  of  pyrimidine  clustered 

with  CH^  and  pyrazine  and  pyrimidine  clustered  with  NH^,  and  benzene 

clustered  with  and  The  pyrimidine-alkane  clusters  are  presented 

7 

for  comparison  with  the  previously  published  pyrazine-alkane  data:  the 
effect  of  the  ring  nitrogen  atoms  on  the  cluster  geometry  can  thereby  be 
evaluated.  The  pyrazine  and  pyrimidine  ammonia  clusters  reveal  the  role  of 
hydrogen  bonding  interactions  in  simple  clusters.  Benzene-water  and-ammonia 
clusters  serve  as  an  example  of  clusters  with  these  more  complicated  solvent 
systems  (i.e.,  two  possible  interaction  potentials)  interacting  with  aromatic 
hydrocarbons.  Theoretical  and  experimental  studies  of  such  a  series  of 
systems  should  eventually  lead  to  a  fuller  understanding  of  solute-solvent 
coordination  structure,  dynamics,  and  the  hydrogen  bonding  interaction. 

Although  extensive  efforts  were  made  and  a  wide  variety  of  experimental 
conditions  explored,  pyrazine  and  pyrimidine  clusters  with  water  were  not 
observed.  Both  FE  and  2-color  70FMS  detection  techniques  were  employed.  A 
broad  feature  (roughly  50  cm-1  FWHMJ  was  observed  in  FE  at  5 HO  cm-1  to  the 
blue  of  pyrazine  0°  transition  for  pyrazine  and  WBter  expanded  with  helium. 

No  signal,  nowever,  was  observed  at  the  pyrazine  water  mass  channel  for 

2-color  T0FN3.  We  must  thus  conclude  that  the  pyrazine  and  pyrimidine  water 

* 

clusters  have  not  been  observed.  Either  the  excited  r.n  state  of  the  cluster 
is  dissociative  or  internal  conversion  or  mtersystec  crossing  is  so  rapid  f : r 
these  clusters  that  the  lifetime  of  the  n*  state  is  greatly  reduced  (£a.  0.1 
ps).  Pynmidine-fluorinated  alcohol  clusters  also  exhibit  broad  features  an: 
reduced  lifetimes. 


.*•  ;• 
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Hydrogen  bonding  interactions  are  known  to  play  an  important  role  in 

intra-  and  mtermolecular  interactions  important  for  secondary  and  tertiary 

molecular  structure,  molecular  dynamics,  and  ionic  and  molecular  solvation. 

Hydrogen  bonding  interactions  have  received  attention  in  previous  supersonic 

molecular  jet  spectroscopic  studies:  these  systems  include  phenol  clustered 

q 

with  various  proton  accepting  molecules,  indole  and  2-ammopyridine  clustered 

with  various  solvents,1^  benzoic  acid  dimers,^1  s-tetrazine  dimers,^ 

13  14 

1 ,4-dihydroxyanthraquinone ,  and  methylsalicylate.  In  most  instances 

blue  shifted  cluster  spectra  are  found  due  to  the  cluster  stabilization  of  tne 

1 5 

ground  state  and  destablization  of  the  excited  state. 

In  this  and  our  other  cluster  studies,  the  experimental  data  are  supple¬ 
mented  by  potential  energy  Lennard-Jones  (LJ)  atom-atom  (b-12-1)  calculations 
of  cluster  structure,  binding  energy,  and  internal  motion.  The  calculations  and 

7 

the  potential  are  discussed  thoroughly  in  our  previous  studies.  The  LJ  calcu¬ 
lations  produce  identical  geometries  but  slightly  lower  binding  energies  than 
the  exp-6  potential  form  used  in  earlier  calculations.  ”  The  LJ  hydrogen 
bonding  (LJ-HB)  potential  form  proves  to  be  more  versatile  than  the  exp-fc  form 
since  many  more  constants  for  different  types  of  atom-atom  interactions  have 
been  independently  reported  for  the  LJ  potential.  ^  None  of  the  potential 
parameters  employed  in  this  work  is  fit  to  the  cluster  data. 

As  is  well  known  and  widely  accepted,  geometries  of  isolated  molecules 
and  clusters  are  best  obtained  through  spectroscopic  observation  of  rotation a! 
structure.  In  fact,  rotational  structure  of  a  number  of  simple  clusters  has 

been  observed  under  molecular  jet  conditions:  s-tetrazine  and  iodine  with  He 
17  18 

and  Ar  and  aniline  with  Ne  and  Ar.  The  resolution  available  to  us  at 
present  is  0.06  cm”1 ;  at  this  resolution  only  rotational  envelopes  are  observ¬ 
able  which  do  not  lend  themselves  readily  to  a  unique  interpretation  of  cluster 
structure.  We  calculate  that  0.005  cm~^  resolution  would  be  required  to 


t 
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resolve  rotational  structure  for  C^N^NH^  etc.,  under  the  restriction 

of  a  rigid  geometry.  We  are  thus  for  the  present  forced  to  employ  less  direct 

methods  to  obtain  cluster  geometry.  Assignment  of  the  spectra  is  accomplished 

through  the  determination  of  ionization  energies,  spectral  shifts,  relative 

intensities,  (molecular)  symmetry  forbidden  cluster  transitions,  and  potential 

calculations.  The  understanding  of  these  more  complex  systems  rests  heavily 

1  -7 

on  the  previous  data  obtained  for  other  clusters.  For  all  systems  dis¬ 
cussed  in  this  paper,  complete  agreement  between  spectroscopic  data,  calculation 
and  results  for  previously  analyzed  solute-solvent  clusters  is  found. 

II.  Experimental  Procedures. 

The  experimental  apparatus  and  procedures  are  similar  to  those  used 
previously  for  the  study  of  vdW  clusters.  The  vacuum  system  consists  of  two 
chambers  witn  a  pulsed  nozzle  and  mass  detection  system  in  the  second  chamber. 

A  skimmer  separates  the  pulsed  nozzle  and  the  time  of  flight  mass  spectrometer. 
The  first  chamber  contains  either  a  pulsed  or  CW  nozzle  the  molecular  beam 
from  which  can  be  taken  into  the  second  chamber  through  a  skimmer.  FE  and  IE 
experiments  are  carried  out  in  the  first  chamber. 

Tne  two  independent  lasers  used  in  the  2-color  TOFMS  experiments  are 
Nd+^/YAG  pulsed  lasers  the  doubled  output  of  which  pumps  two  dye  lasers. 

The  dye  laser  output  can  be  nixed  with  the  1.064  pm  Nd^/YAG  fundamental, 
frequency  mixed  and  doubled,  or  just  doubled  using  various  nonlinear  KIP 
crystals.  The  laser  output  can  be  extended  from  greater  than  4.5  yr  to 

-0.215  um.  Ore  laser  is  employed  to  excite  the  cluster  to  its  first  excited 

•  • 

nir  or  wit  electronic  state  and  the  second  laser  then  ionizes  this  cluster 
starting  from  the  £ ^  vibronic  manifold.  The  maximum  ionization  energy 
achievable  with  this  second  laser  is  roughly  46,500  cm"1. 

The  solute  or  cluster  chronophore  is  typically  placed  in  an  inline  trap 


or  filter  cup  directly  behind  the  pulsed  valve.  Water  is  placed  in  a  trap 


before  the  valve  and  solute;  the  helium  carrier  gas  passes  over  both  materials 
and  into  the  valve.  Gaseous  solvents  are  premixed  with  the  helium  carrier  gas 
in  a  holding  tank  (-2000  psi)  at  concentrations  varying  from  2.0  to  0.1  mole/ 
mole  percent. 

The  LJ  potential  function  (6-12-1)  with  the  additional  HB  form  is 

7 

described  in  detail  in  a  previous  publication.  Table  I  contains  a  list  of 
the  previously  unreported  constants  employed  with  this  potential  form.  Pyrazme 
and  pyrimidine  structures  used  in  the  calculations  are  obtained  from  reference 
19. 

III.  Results. 

This  section  contains  the  experimental  and  calculational  results  for 
the  various  clusters  investigated.  We  first  present  pyrimidine  with  methane 
and  ethane  for  comparison  with  the  previously  reported  pyrazine  clusters. 1 
Pyrimidine  and  pyrazine  ammonia  clusters  are  then  discussed  and  benzene  water 
and  ammonia  clusters  are  presented  for  a  comparison  with  the  N-heterocyclic 

1  -n 

systems.  Based  on  previous  experience  with  a  number  of  different  clusters, 
we  have  not  made  an  extensive  experimental  study  of  the  binding  energies  of 
these  clusters.  We  rely  on  the  calculations  which  have  always  fallen  within 
tne  range  bracketed  by  the  experimental  2-color  TOFKS  deta.  Preliminary 
checks  made  on  both  pyrazine  and  pyrimidine  clusters  are  in  agreement  with  the 
calculations.  Vibrational  vdW  modes  observed  in  the  spectra  will  not  be 
assigned  in  this  publication.  Future  publications  will  assign  then  as  bends, 
stretches,  torsions  and  combinations  utilizing  a  normal  coordinate  calcula- 
tional  analysis. 

A.  Pyrimidine  -  methane. 

Figure  1  and  Table  II  present  the  data  for  the  pyrimidine  (CK  ) 
and  (CW^)^  complexes  taken  near  the  pyrimidine  0°  transition  (31073.0 
cm  1).  The  cluster  of  pyrimidine  (CH^)^  has  a  spectral  shift  of  -5^.6 
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cm”1  and  a  low  frequency  mode  at  4.5  cm”1  from  this  origin.  These  spectra 

are  obtained  by  2-color  TOFMS.  Hints  of  other  vdW  modes  can  also  be  seer  in 

the  trace  in  figure  1  but  we  are  hesitant  to  report  such  weak  transitions. 

Pyrimidine  (CH^)^  spectra  clearly  show  two  clusters,  similar  to  previously 

1-7  -1 

reported  clusters  for  other  aromatic  systems.  The  feature  at  -112.1  cm 

in  figure  1  is  assigned  as  the  0°  transition  of  the  isotropic  (symmetric), 

additive  shift  cluster  and  the  feature  at  -47.2  cm”1  is  attributed  to  the  0° 

o 

transition  of  the  anisotropic  (asymmetric)  cluster  with  both  methanes  on  the 
same  side  of  the  aromatic  ring.  Note  that  without  both  mass  and  energy  resolu¬ 
tion,  the  spectra  of  pyrimidine  (CH^  and  ( CH^ )2  would  not  be  resolved  and 
the  clusters  could  not  be  separated  and  uniquely  identified.  A  vdW  mode  at  5.1 
cm”1  from  the  isotropic  cluster  origin  is  observed. 

Potential  energy  calculations  using  LJ  potentials  for  these  clusters 
generate  geometries  and  binding  energies  comparable  to  those  previously  resorted 
for  other  aromatic-alkane  systems  and  in  complete  accord  with  the  above  experi¬ 
mental  findings  (see  fig.  2).  Calculations  for  pyrimidine  (CH^  clusters 
yield  only  one  geometry  for  which  the  methane  is  coordinated  with  the  aromatic 

w-system  of  the  pyrimidine  ring.  The  calculated  binding  energy  for  this  cluster 

-1  0 
is  514  cm  .  The  methane  carbon  atom  is  above  the  ring  at  3*5  A  and  is  shifted 

-0.1  X  from  the  ring  center  toward  the  nitrogen  atoms.  The  three  hydrogen  atoms 

o 

of  methane  that  point  down  toward  the  ring  are  equidistant  from  the  ring  at  3.1 
two  of  these  hydrogens  point  directly  at  the  ring  nitrogen  atoms.  Again  in 
agreement  with  the  main  experimental  observations  for  these  clusters,  the  calcu¬ 
lations  for  pyrimidine  (CH^)^  clusters  yield  two  distinct  geometrical  arrange¬ 
ments.  The  isotropic  cluster  has  a  calculated  binding  energy  of  1029  cm”1  and 
the  anisotropic  cluster  has  a  calculated  binding  energy  of  879  cm”1. 
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B.  Pyrimidine  -  Ethane. 

The  pyrimidine  spectrum  is  quite  complicated,  consisting 

of  a  number  of  low  intensity  features  and  an  intense  feature  at  -60.7  cm-1 

with  respect  to  the  pyrimidine  origin  (see  figure  3)»  We  might  expect  that 

little  vdW  vibronic  intensity  would  be  observed,  based  on  pyrazine  and  other 

pyrimidine  clusters.  In  order  to  begin  to  interpret  these  features  we  must 

21 

consult  the  vapor  phase  room  temperature  pyrimidine  monomer  spectrum.  In 
these  reports,  features  at  -156  cm  1  and  +22  cm  1  are  assigned  as  the  16a| 
and  16bj  sequence  bands,  respectively.  The  feature  that  appears  in  figure  5 
at  —39 • 1  cm  1  may  be  assigned  as  the  16b1  sequence  band  of  the  pyrimidine 
^2^6^1  c^u8l:er  built  on  the  intense  -60.7  cm*1  cluster  origin.  If  this 
identification  is  correct,  then  the  clusters  of  pyrimidine  are  6ot 

(Tyit  »  260  K)  and  the  features  in  figure  3  and  table  II  at  -153*5,  -99*3, 

-86.9,  and  -75«4  cm  1  may  well  be  hot  bands  associated  with  the  cluster  origins 
at  -71.4,  -60-7  and  -52.7  cm-1.  Apparently  the  cluster  formation  process  for 
pyrimidine  ethane  tends  to  warm  the  cluster.  Table  III  summarizes  these 
results. 

Supersonic  expansion  of  pyrimidine  apparently  does  not  produce  the  expected 
cooling  for  vibrational  modes  1 6b ^  and  16a^.  Vibrational  temperatures  for  the 

n  i  j 

16b^  mode  have  been  reported  to  be  in  excess  of  200  K.  Ito  and  coworkers 
have  also  observed  the  16b|  transition  for  pyrimidine  clustered  with  argon  and 
nitrogen.  In  the  present  work,  changes  in  backing  pressure  from  10  to  120  psi 
do  not  change  the  relative  band  intensities:  apparently  these  modes  present  a 
bottleneck  for  vibrational  cooling. 

As  can  be  seen  in  figure  4,  three  different  configurations  are  calculated 
to  be  stable  for  the  pyrimidine  (C2Hg)1  cluster.  Configuration  I  has  the 
long  axis  of  ethane  perpendicular  to  the  plane  of  the  pyrimidine  ring. 
Configurations  II  and  III  have  the  ethane  molecule  long  axis  more  or  less 


parallel  to  the  plane  of  the  ring;  for  configuration  II,  the  axis  of  ethane 
lies  between  a  nitrogen  and  a  carbon  and  for  configuration  III  this  axis  lies 
between  two  carbon  atoms.  In  each  instance  a  CH^-group  lies  more  or  less  over 
the  ring  center. 

C.  Pyrimidine  -  Ammonia. 

The  pyrimidine  (NH^)^  2-color  TOFMS  spectrum  in  the  range  300  to 
500  cm-1  to  the  blue  of  the  pyrimidine  origin  is  presented  in  figure  5.  The 
three  features  are  associated  with  the  0°  transitions  of  the  pyrimidine 
(NH^)1  clusters.  The  lack  of  significant  vdW  vibrational  mode  intensity 
indicates  that  the  ground  and  excited  state  vdfo  potentials  are  nearly 
identical.  The  large  cluster  blue  shift  implies  a  strong  hydrogen  bonding 
interaction  between  the  pyrimidine  and  the  ammonia.  The  spectrum  of  figure  5 
is  taken  with  an  ionization  energy  of  45,110  cm-1;  lowering  the  ionization 
energy  to  ca.  42,000  cm-1  the  two  nearly  degenerate  features  at  +367 
cm-1  decrease  in  intensity  much  more  rapidly  than  the  feature  at  +496 
cm-1 .  Table  IV  gives  the  energies  and  assignments  for  the  pyrimidine 
(NK^)^  clusters.  Based  on  the  idea  that  the  three  features  in  this  spectrum 
are  associated  with  three  different  pyrimidine  (NH^)1  cluster  geometries, 
the  two  nearly  isoenergetic  configurations  must  be  quite  similar.  The  feature 
at  +496  cm-1  must  correspond  to  a  different  (more  hydrogen  bonded)  geometry. 

Three  different  configurations  are  calculated  for  the  pyrimidine  (NH^)1 
cluster  using  the  LJ-HB  potential  form  (see  Table  I  and  Reference  7).  Two  of 
these  configurations  have  the  ammonia  above  the  pyrimidine  plane  and  the  third 
configuration  has  the  NH^  molecule  not  above  the  ring.  These  clusters  are 

presented  m  figure  6  along  with  the  calculated  cluster  binding  energies.  Con- 

o 

figuration  I  has  the  ammonia  R-atom  3*2  A  above  the  pyrimidine  ring  with  the 

three  ammonia  hydrogen  atoms  pointing  down  toward  the  ir -cloud.  The  two  closest 

o 

hydrogen  atoms  of  the  ammonia  in  this  configuration  are  2.7  A  from  the  pynmi- 


dine  plane.  In  configuration  II,  the  ammonia  N-atom  is  3*2  X  above  the  ring 

plane  and  points  in  the  general  direction  of  the  ring  N-atoms.  The  closest  two 

o 

ammonia  H-atoms  to  the  ring  are  at  2.7  A.  This  cluster  has  a  somewhat  lower 
binding  energy  than  cluster  I  (667  vs  689  cm-1).  In  the  third  pyrimidine 
(NHj)^  configuration,  the  ammonia  molecule  does  not  reside  over  the  pyrimidine 

ring:  the  NH^  has  a  N-H  bond  in  the  plane  of  the  pyrimidine  ring  with  the  H- 

o 

atom  pointing  to  one  of  the  pyrimidine  nitrogens,  2.23  A  from  it.  The  ammonia 
nitrogen  is  displaced  in  the  xy  plane  by  3«3  $  and  2.8  X  along  the  x-axis  and 
y-axis  respectively,  as  shown  in  figure  6,  configuration  III.  This  configura¬ 
tion  has  a  calculated  binding  energy  of  537  cm-1 . 

D.  Pyrazine  -  Ammonia. 

The  pyrazine  (NH_),  spectra  in  the  0°,  10a1  and  6a1 

J  1  0  0  0 

regions  are  presented  in  figure  7-  One  first  notices  the  significant  dif¬ 
ference  between  these  data  and  those  of  pyrimidine  (NH^.  The  spectra 
all  strongly  suggest  that  only  one  configuration  is  present  for  the  pyrazine 
(NH^)1  system.  The  intense  feature  in  the  0°  spectrum  is  the  cluster 
origin  at  +117  cm-1  from  the  pyrazine  0°  transition.  The  remaining 

features  of  this  spectrum  are  vdW  vibrational  modes  of  the  cluster.  They 

20 

will  be  analyzed  in  a  future  publication.  Table  V  gives  the  energies  and 

the  features  observed  in  the  spectrum.  In  this  case,  the  ground  and  excited 

state  potential  surfaces  must  be  significantly  different. 

The  cluster  10a1  transition  is  quite  different  from  the  cluster  0°  or 
o  ’  o 

6a1  transitions.  From  this  one  concludes  that  strong  vdW-internal  mode 

1  22 

coupling  exists  for  the  10a  out  of  plane  ring  mode.  The  vdW  overtones 

and  combination  bands  extend  to  more  than  180  cm”1  from  the  10a1  cluster 

o 

origin  feature.  This  is  particularly  striking  in  comparison  with  the  6a1 


vibrornc  band. 


Only  one  configuration  is  calculated  for  the  pyrazine  (NH^)^  cluster 
using  the  LJ  plus  HB  potential  function,  in  agreement  with  expectations  from 
the  spectra.  The  ammonia  molecule  hydrogen  bonds  to  the  ring  nitrogens, 
through  two  ammonia  hydrogens.  The  nitrogen  atom  of  the  ammonis  molecule  is 
3.2  $  above  the  ring  plane  with  all  three  hydrogens  pointing  toward  the  ring. 
This  geometry  is  depicted  in  figure  8.  The  calculated  binding  energy  of  the 
cluster  is  677  cm-1  . 

E.  Benzene  -  Ammonia. 

The  origin  and  6^  transitions  of  the  CgHg(NH^)^  cluster  are 

observed  in  2-color  TOFMS.  One  can  immediately  conclude  that  at  least  one 

configuration  of  this  cluster  does  not  retain  the  benzene  3-fold  axis.  The 

spectra  are  traced  in  figure  9-  The  spectra  for  both  transitions  are  far  more 

complicated  than  any  other  cluster  previously  reported  from  our  laboratory. 

Regions  with  similar  structure  can  be  found  in  the  0°  and  61  spectra;  in 

o  o 

particular,  the  regions  around  +20  cm-1 ,  -20  cm-1  and  -60  cm-1  in  each 

spectra  bear  some  resemblance  to  one  another.  Nonetheless,  features  appearing 

in  one  spectrum  do  not  appear  in  the  other,  and  therefore  at  least  two 

clusters  of  different  geometry  are  probably  responsible  for  the  61 

o 

transition.  The  0°  transition  could  arise  from  one  reduced  symmetry  cluster. 
Since  very  little  £  priori  spectroscopic  analysis  seems  possible 
in  this  situation  at  the  present  time,  great  stocF  must  be  placed  in  the 
calculations. 

Two  different  configurations  are  calculated  for  the  C,H,(NH,), 

bo  3  ' 

cluster  and  both  of  them  have  the  NH^  molecule  placed  over  the  ring. 
Configuration  I  in  figure  10  preserves  the  benzene  axis  and  configuration 
II  does  not.  The  high  symmetry  cluster  has  the  N-atom  3*3  X  from  the  ring 
plane  and  the  three  H-atoms  of  NH^  are  2.9  $  above  the  ring  plane.  Configura¬ 
tion  II  has  the  N-aton  3*3  X  above  the  plane  and  slightly  shifted  from  the  ring 


center.  Two  H-atoms  of  the  NH^  point  towards  C-C  bonds  and  one  points  away 
from  the  ring  plane.  The  NH^  H-atoms  in  this  configuration  are  at  2.9  X  (two 
of  them)  and  4.3  X  from  the  ring.  The  calculated  binding  energies  of  these  two 
configurations  are  given  in  figure  10.  Configuration  II  alone  must  generate 
-.he  0°  spectrum  in  figure  9« 

F.  Benzene  -  Water. 

The  C6H6(H20)1  cluster  0°  and  6^  spectra  are  presented  in 

figure  11.  They  are  strikingly  different  from  those  of  C^H^(NH^)1 

(figure  9)  but  bear  a  strong  resemblance  to  the  pyrazine  (NH^)^  spectra 

displayed  in  figure  7,  as  well  as  spectra  of  other  systems  studied  in  our 
1  -7 

laboratory.  In  particular,  a  very  clear  vibronic  progression  is  present 

for  the  0°  spectrum  and  will  be  analyzed  in  a  future  publication. ^  Table  VI 

gives  these  features  and  their  energies.  The  61  spectrum  of  CcH, (H„0),  has 

0  o  b  d.  i 

a  different  shift  (+50  vs  +85  cm-1)  and  a  different  general  pattern  than 
the  0°.  In  addition,  the  6^  origin  is  split  by  roughly  1  cm*1 .  Consider¬ 
able  vibronic  coupling  must  occur  between  the  in  plane  carbon-carbon  stretch 
61  and  the  vdW  modes.  The  CgHg(H20)^  cluster  2-color  TOFMS  spectra  are 

observed  only  at  a  higher  ionization  energy  than  required  to  observe  the  61 

o 

spectrum  of  bare  C^H^.  The  C^H^(H20)^  spectra  depicted  in  figure  11  are 
taken  with  an  ionization  laser  energy  of  44,480  cm-1  compared  to  36,100  cm-1 
for  the  benzene  monomer.  The  ionization  energy  for  the  benzene  water  cluster 
is  3300  cm-1  higher  than  that  of  the  benzene  monomer. 

Only  one  geometry  is  calculated  for  the  C^H^(H20)1  cluster  employing 
the  LJ  potential.  The  configuration  and  binding  energy  are  given  in  figure 
12.  The  oxygen  atom  of  water  is  roughly  centered  over  the  ring  at  3*2  A  above 
it  and  the  water  hydrogen  atoms  are  at  3.0  X  above  the  ring.  The  calculated 
water  benzene  binding  energy  is  505  cm”1 . 


II 
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IV.  Discussion 

Cluster  geometry  is  determined  through  analysis  of  2-color  TOFMS  data  for 
individual  clusters  and  through  calculations  of  cluster  geometry  and  binding 
energies  using  an  augmented  LJ  potential.  Calculated  binding  energies  have 
always  been  roughly  bracketed  by  the  experimentally  observable  range  defined 
by  two  solute  molecule  vibrations  (e.g.t  520  to  850  cm  for  the  1:1  complex). 

Before  discussing  the  details  of  the  clusters  observed,  some  general 
remarks  are  in  order.  First,  one  should  take  note  of  the  great  apparent 
differences  between  the  spectra  of  the  various  clusters  studied:  pyrazine  and 
pyrimidine  hydrocarbon  spectra  show  little  vdW  vibronic  structure  but  intense 
origins;  the  pyrazine  (NH^)^  spectrum  displays  elaborate  and  well  developed 
vibronic  progressions  with  a  +117  cm-1  cluster  shift  for  a  single  cluster; 
the  pyrimidine  (NH^)^  spectra  consist  only  of  intense  origins  for  three 
large  shift  (+366,  +368,  +496  cm-1)  configurations  with  no  vdW  vibronic 
development;  the  benzene  (NH^)1  spectra  are  red  shifted,  and  too  complex  to 
interpret  without  further  calculations;  and  the  benzene  (H20)1  spectra  are 
blue  shifted  with  extensive  vdW  structure.  Second,  and  perhaps  even  more 
astonishing,  the  LJ-HB  (where  appropriate)  potential  calculations  parallel  and 
reinforce  these  differences  in  all  cases.  That  is,  for  example,  the  calcula¬ 
tions  suggest  one  hydrogen  bonded  configuration  for  pyrazine  (NH^)1  but  three 
hydrogen  bonded  configurations  for  pyrimidine  (Nh^,  in  agreement  with  the 
straight-forward  interpretation  of  the  spectra.  Third,  the  binding  energies  of 
the  clusters  seem  relatively  insensitive  to  the  detailed  configuration  of  the 
cluster.  Fourth,  spectral  shifts  are  found  to  be  a  sensitive  function  of  the 
detailed  geometry  of  the  cluster.  Proximity  to  the  n-system  is  important  for 
cluster  red  shifts,  while  hydrogen  bonding  yields  in  general  cluster  blue  shifts 
with  respect  to  the  solute  monomer  origin. 


A.  Pyrimidine  -  Methane. 

In  the  pyrimidine  (CH^)^  cluster  the  CH^  molecule  is  situated 

above  the  pyrimidine  ring  coordinated  to  the  it-cloud  of  the  aromatic  ring. 

The  cluster  has  a  simple  spectrum  with  little  vdW  vibronic  intensity  following 

the  pyrimidine  vibronic  origins.  The  cluster  spectrum  is  red  shifted, 

indicating  that  the  excited  state  cluster  is  more  tightly  bound  than  the 

ground  state  cluster  by  roughly  60  cm”^.  The  overall  appearance  of  the 
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spectrum  is  similar  to  that  of  pyrazine  methane. 

The  additive  shift  features  in  the  spectrum  of  pyrimidine  (CH^)^  are 
attributed  to  the  isotropic  (symmetrical)  geometry  with  a  methane  molecule  on 
either  side  of  the  pyrimidine  ring.  The  feature  at  -47.2  cm-1  in  the 
pyrimidine  (CH^)^  spectrum  is  attributed  to  the  anisotropic  configuration 
with  both  CK^  molecules  on  the  same  side  of  the  ring.  In  this  asymmetric 
geometry,  one  methane  is  more  or  less  above  the  ring,  and  responsible  for  most 
of  the  CH^-ir-cloud  interactions,  and  the  other  methane  molecule  is  off  the 
ring  interacting  primarily  with  the  first  methane,  contributing  little  to  the 
cluster  spectral  shift. 

The  overall  behavior  of  this  cluster  system  with  regard  to  geometry, 

binding  energy,  and  cluster  population  in  the  beam  is  very  similar  to  that 

1  -7 

found  for  benzene,  toluene  and  pyrazine  methane  species. 

B.  Pyrimidine  -  Ethane. 

The  pyrimidine  cluster  is  similar  to  the  pyrazine 

n 

cluster:  both  clusters  have  three  geometrical  configurations  and 

C  O  1 

each  configuration  of  the  two  clusters  has  a  similar  binding  energy.  These 
configurations  are  also  similar  to  those  of  the  benzene  clU8t@r 

with  the  exception  that,  in  the  N-heterocyclic  systems,  two  parallel 
orientations  of  the  long  axis  with  respect  to  the  ring  plane  are  now 

possible. 

13 


The  identification  of  calculated  geometries  with  the  three  distinct 
spectroscopic  features  is  of  course  tentative  but  can  be  pursued  in  the  spirit 

of  the  arguments  and  correlations  employed  with  the  benzene  and  pyrazine 

1-7  -1 

systems.  Referring  to  figures  3  and  4,  the  feature  at  -71.4  cm  can  be 

associated  with  configuration  1,  the  intense  feature  at  -60.7  cm”1  can  be 
associated  with  configuration  II,  and  the  -52.7  cm”1  feature  is  associated 
with  configuration  III.  The  lack  of  vdW  vibronic  structure  for  these  transi¬ 
tions  must  be  due  to  the  similarity  between  the  ground  and  excited  state 
potentials  for  the  clusters. 

C.  Pyrimidine  -  Ammonia. 

The  pyrimidine  (NH^)^  clusters  have  a  unique  spectrum  which  can 
only  be  interpreted  as  due  to  three  distinct  configurations  with  no  vdW  mode 
progression  intensity  following  the  well  defined  origins.  Calculations,  as 
pointed  out  previously,  give  exactly  these  conclusions  and  identification  of 
origins  in  the  spectrum  with  configurations  seems  straightforward.  Configur¬ 
ations  I  and  II  of  figure  6  are  associated  with  the  features  at  366  and  368  cm” 
in  figure  5*  These  two  configurations  are  quite  similar  and  have  less  hydrogen 
bonding  interaction  than  the  more  blue  shifted  single  feature  at  496  cm”1 . 

The  large  spectral  blue  shifts  of  -365  and  -500  cm”1  must  arise  from  the  strong 
hydrogen  bonding  interactions.  While  none  of  these  observations  seem  par¬ 
ticularly  striking  in  and  of  itself,  in  comparison  with  the  pyrazine  and 

benzene  (NH^)1  results,  they  are  surprising;  these  will  be  discussed  below. 

B.  Pyrazine  -  Ammonia. 

The  pyrazine  (NH^)^  cluster  spectrum  is  completely  different 

from  that  of  pyrimidine  (NH, ),.  The  cluster  0°  is  shifted  +117  cm”1 

3  '  o 

from  the  pyrazine  0®,  only  one  cluster  geometry  is  present,  and  a  rather 

extensive  vdW  vibronic  structure  is  built  upon  the  origin.  The  6a1  spectrum 

0 

is  quite  similar.  Two  intense  vdW  vibronic  transitions  are  associated  with 
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these  cluster  transitions.  The  10a^  vibration  region  of  the  pyrazine  (NH^)1 
cluster,  on  the  other  hand,  looks  quite  different  from  these  other  features  and 
vibronic  interactions  between  the  out  of  plane  10a^  carbon-carbon  mode  and  the 
vdW  modes  are  quite  evident  in  the  overall  vibronic  intensity  pattern  in  the 
10a^  region. 

Calculations  predict  only  one  configuration  for  this  cluster  system.  The 
fact  that  the  potential  energy  calculations  can  accurately  parallel  the 
spectral  data  for  pyrimidine  and  pyrazine  ammonia  clusters  which  are  clearly 
so  different,  gives  us  a  high  degree  of  confidence  in  the  calculational 
process,  the  binding  energies,  and  the  potential  form  accuracy. 

E.  Benzene  -  Ammonia. 

Assignment  of  the  benzene  (NH^)^  spectra  has  not  been  attempted 
as  yet  because  the  spectra  are  too  complicated.  We  present  them  only  as  part 


of  the  general  picture  indicating  what  spectra  of  relatively  simple  clusters 
of  such  systems  can  be  like.  The  benzene  (NH^)^  clusters  yield  much  more 
complicated  spectra  than  either  pyrazine  or  pyrimidine  (NH^)^  clusters  do. 

In  the  benzene  (NH^)^  cluster,  strong  vdW  vibronic  interactions  must  be 
important  for  the  cluster  transition  intensity. 


Both  calculated  configurations  (figure  10)  contribute  to  the  6 


1 


transition  but  only  configuration  II  generates  the  0  spectrum.  Note  too, 

o 


that  both  configurations  most  likely  generate  a  red  shift. 

F.  Benzene  -  Water. 

Any  CgHgtH^O)^  cluster  will  in  principle  generate  a  0^ 
transition.  The  0°  spectrum  of  (figure  11)  thus  strongly 

suggests  that  only  one  cluster  geometry  is  realized  for  this  system.  The 


)° 

o 


0°  transition  is  much  like  the  pyrazine  (NH^)1  spectrum. 


Well  developed  vdW  vibronic  features  are  observed.  The  6^  spectrum  implies 


strong  vibronic  mixing  between  the  in  plane  carbon-carbon  deformation  6  and 
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the  vdW  bends  and  torsions:  the  cluster  shift  is  different  for  61 .  as  is 

o 

the  intensity  pattern.  The  cluster  shift  at  the  0°  transition  is  +85  cm-1 

which  indicates  that  the  cluster  excited  state  is  destablized  with  respect  to 

the  ground  state.  The  blue  shift  may  be  related  to  the  unique  hydrogen 

bonding  capabilities  of  the  ^0  molecule  with  the  ft-system  of  benzene.  The 

calculated  binding  energy  is  probably  -50  cm”1  low  for  this  cluster  because 

the  6^  transition  at  0°  plus  520  cm”1  is  observed.  We  have  previously 

noted  that  the  LJ  potential  binding  energy  is  roughly  50  cm”1  low  compared 

7 

to  exp-6  and  experimental  values. 

An  infrared  study  of  the  1:1  benzene-water  complex  in  an  argon  matrix  has 
23 

been  reported.  The  work  suggests  that  the  water  molecule  hydrogen  bonds 
to  the  benzene  ir-system  in  a  manner  nearly  identical  to  that  found  in  the 
calculations  presented  in  figure  12. 

V.  Conclusions. 

The  clusters  studied  in  this  work  fall  into  two  broad  categories:  a 
conventional  set  containing  pyrimidine  hydrocarbon  clusters,  the  spectra  of 
which  are  quite  similar  to  those  of  other  aromatic  and  pyrazine  hydrocarbon 
systems,  and  benzene,  pyrimidine,  and  pyrazine  ammonia  and  benzene  water 
clusters,  the  spectra  of  which  are  all  unique  and  surprisingly  erratic.  In 
the  latter  grouping,  spectral  cluster  shifts  range  from  -100  to  +500  cm”1, 
vdW  vibronic  spectra  range  from  nonexistent  to  intense,  vdW  modes  can  be  highly 
perturbing  to  the  solute  vibronic  structure  and  energy,  and  the  nunber  of 
cluster  configurations  varies  from  one  to  three  in  an  apparently  arbitrary 
fashion.  These  differences  notwithstanding,  the  Lennard-Jones  (6-12-1)  poten¬ 
tial,  augmented  appropriately  with  hydrogen  bonding  interactions  (10-12)  as 
required,  always  gives  geometry  and  binding  energy  results  that  are  in  complete 
agreement  with  the  spectra  as  far  as  the  comparison  can  be  made  (i.e.,  symmetry, 
numbers  of  configurations,  red  and  blue  shifts  with  regard  to  hydrogen  bonding, 


etc).  The  atom-atom  LJ  potential  form  has  been  chosen  for  these  calculations 
because  a  large  number  of  parameters  for  different  types  of  atoms  are  available 
in  the  literature. 

The  benzene  (NH^)^  and  (l^O)^  spectra  are  quite  different  from  one 
another.  The  shifts  for  these  two  clusters  and  their  geometries  seem  to  empha¬ 
size  the  importance  of  hydrogen  bonding  in  the  benzene  (H^O)^  cluster. 

Pyrazine  and  pyrimidine  water  clusters  are  not  found  in  these  studies 
although  they  have  been  extensively  investigated.  These  clusters  are  not 
observed  perhaps  because  their  excited  states  are  dissociative,  but  more 
likely  because  of  rapid  excited  state  irtersystem  crossing  and/or  internal 


conversion. 
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Table  II 


Observed  peaks  in  the  spectra  of  pyrimidine-methane  clusters. 


Species 

Energy. 1 
(vac.  cm  ) 

Energy  Relative 
to  Pyrimidine  0o 
(cm“l) 

Energy  Relative 
to  Cluster  0° 
(cm-1 ) 

Assignment 

C4H4N2(CVl 

31016.4 

-56.6 

0 

o  o 
o 

31020.9 

-52.1 

4.5 

C4H4N2(CH4}2 

30960.9 

-112.1 

0 

iso 

o  o 
o 

30966.0 

-107.0 

5.1 

31025.8 

-47.2 

0 

aniso 

0° 

o 


a)  The  -112*1  cm-1  shift  is  associated  with  two  methanes  added  symmetrically  above  and 
below  the  pyrimidine  ring  (isotropic),  as  shown  in  Figure  2.  The  -47.2  cm-1  shift 
is  associated  with  the  anisotropic  configuration,  as  shown  in  Figure  2. 
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Observed  peaks  in  the  spectra  of  pyrimidine-ethane  clusters. 
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Table  V 


Pyrazine 


Energy  1 
(vac.  cm”  ) 

Energy  Relative 
to  corresponding  1 

pyrazine  feature  (cm  ) 
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Benzene 


Energy 
(vac.  cm*  ) 


Energy  Relative 
to  corresponding  1 
benzene  feature  (cm*  ) 
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water  "feature  (cm*  ) 
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FIGURE  CAPTIONS 


Figure  1  Two-color  time  of  flight  mass  spectra  (2-color  TOFMS)  of 

pyrimidine  (CH^)1  and  pyrimidine  (CH^)^  in  the  region 
of  the  pyrimidine  origin  (31073.0  cm-1). 

Figure  2  Minimum  energy  configurations  and  binding  energies  for 

pyrimidine  (CH^)^  and  pyrimidine  (CH^^  as  obtained  using 
the  LJ  potential  calculation  described  in  the  text. 

Figure  3  Two-color  TOFMS  of  pyrimidine  ( ^  in  the  region  of  the 

pyrimidine  origin. 

Figure  4  Minimum  energy  configurations  and  binding  energies  for 

pyrimidine  ( ^  as  obtained  using  a  LJ  potential 
calculation. 

Figure  5  Two-color  TOFMS  of  pyrimidine  (NH^)1  in  the  region  300-500 

cm-1  to  the  blue  of  the  pyrimidine  origin.  An  insert  is 
shown  of  the  first  two  features  with  an  expanded  scale. 

Figure  6  Minimum  energy  configurations  and  binding  energies  for 

pyrimidine  (NH^)^  as  obtained  using  a  LJ  plus  HB  potential 
calculation. 

Figure  7  Two-color  TOFMS  of  pyrazine  (NH^)^  at  the  origin,  10a^  and 

6a^  regions  are  shown.  The  cluster  origin  (0  cm  ^ )  is 
♦117  cm-1  to  the  blue  of  the  pyrazine  origin. 


Figure  8 


Minimum  energy  configuration  and  binding  energy  for  pyrazine 


Figure  9 


Figure  10 


Figure  1 1 


Figure  12 


(NH^)^  as  obtained  using  a  LJ  HB  form  potential  calculation. 

Two-color  TOFMS  of  benzene  (NH^)^  in  the  region  of  the  benzene 

0°  and  6^  transitions, 
o  o 

Minimum  energy  configurations  and  binding  energies  for  benzene 
(NH^)1  as  obtained  using  a  LJ  potential  calculation. 

Two-color  TOFMS  of  benzene  in  the  region  of  benzene 

origin  0°  and  61 .  An  expanded  scale  insert  is  shown  of 
o  o 

the  6^  origin  of  the  cluster. 

Minimum  energy  configuration  and  binding  energy  for  benzene 
(h2o)  1  as  obtained  using  a  LJ  potential  calculation. 
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